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Abstract. Competition between the vortex lattice and a lattice of asymmetric 
artificial defects is shown to play a crucial role in ratchet experiments in 
superconducting films. We present a novel and collective mechanism for current 
reversal based on a reconfiguration of the vortex lattice. In contrast to previous 
models of vortex current reversal, the mechanism is based on the global response 
of the vortex lattice to external forces. 
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1. Introduction 

The term ratchet effect refers to a net flow of particles induced by asymmetric 
potentials, in the absence of non-zero average forces. This effect has drawn the 
attention of many researchers from very different fields [Tj. It is believed to play 
a role in protein motors [2] and it has been used to design synthetic molecular motors 
[21 2] and control the motion of colloidal particles [5] . 

In the context of superconductivity, vortices in superconducting films [51 [7| and 
Josephson junctions [HIS] were soon proposed as implementations of the ratchet effect. 
The first experimental evidence of rectification of vortices was reported by Villegas et 
al. [10] . In this work, a Nb superconducting film was grown on an array of triangular 
Ni nanoislands, acting as pinning sites for the vortices. When the vortex lattice is 
driven by an input ac-current, the asymmetric geometry of the pinning potential leads 
to a net flow of vortices, resulting in an output dc-voltage. After this seminal work, 
considerable theoretical [TTl [T^l IS] and experimental [2] research has been done in 
order to clarify ratchet mechanisms or use them to control the motion of vortices in 
fllms and Josephson arrays. 

An interesting feature of this experimental setup is the appearance of current 
reversals, i.e., the rectiflcation polarity changes for certain values of the applied 
magnetic field and current strength. This behavior has been explained in Ref. [10] 
by means of a one dimensional model, where current reversal is due to the interplay 
between vortices pinned in the asymmetric centers (triangles) and interstitial vortices 
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placed among the triangles. However, in that model the interaction between pinned 
and interstitial vortices was replaced by an effective ratchet potential affecting only 
to interstitial vortices. Beside this effective potential, both pinned and interstitial 
vortices were treated as independent particles. 

In the present paper, we review the experimental results of Ref. [10] of rectification 
and current reversal using magnetic defects, and present experimental results showing 
that the current reversal persists with non-magnetic defects. Furthermore, we 
introduce a new theoretical approach to the current reversal considering interaction 
between vortices. Furthermore, we introduce a new theoretical approach to the current 
reversal considering interaction between vortices. In the absence of pinning potentials, 
it is well known that the ground state induced by vortex interaction is a triangular 
Abrikosov lattice [15]. Our aim is to analyze how this vortex lattice is distorted 
by the lattice of asymmetric defects. The competition between the two lattices will 
be relevant when the magnetic field creates a few vortices per pinning site and for 
long penetration depth. A, which is the case for temperatures close to the critical 
temperature, as in the experiments reported in [10] . Furthermore, we will show that 
this vortex lattice distortion can induce current reversal but nevertheless the lattice 
exhibits a net motion, i.e., both pinned and interstitial vortices move in the same 
direction. Therefore, our picture of vortex rectification and current reversal differs 
both from the one presented in Ref. [TD], where the latter is induced only by the 
motion of interstitial vortices, and that in Refs. [161 [H] which deal with systems 
of particles where adding particles of an auxiliary species, that interact with both 
the primary particles of interest and the substrates, the particles could move either 
together or in opposite directions depending upon the strength of the interaction, and 
whether it is attractive or repulsive. 

Moreover, in our new approach, the superconductor ratchet becomes a collective 
ratchet. Collective effects in ratchets have been shown to yield new interesting 
phenomena such as instabilities and hysteresis [2], absolute negative resistance [17] 
and inefficiency of some feedback control protocols [H]. Moreover, many biological 
motors, such as kinesins and myosins work in a collective way [5]. Collective effects 
have been also considered in the context of superconducting vortices by Reichhardt 
et al. [19] and de Souza Silva et al. [20]. Also Shalom and Pastoriza suggested that 
collective effects are responsible for current reversal in Josephson junction arrays |21j . 

In our approach we show that the ratchet reversal corresponds to a reconfiguration 
of the vortex lattice at different fields and drives. Similar arguments about the 
importance of the vortex- vortex interactions in current reversal mechanisms were made 
in Ref. [19] . In their simulations they observe how interstitial vortices may help depin 
vortices at the tip of the triangular defect due to thermal fiuctuations but not the 
ones at the base of the triangle, leading to a negative current at low drives. Also, the 
current reversal is explained by vortex-vortex interactions in the framework of a one- 
dimensional model in [20j . However, none of these works discuss the aforementioned 
competition between the Abrikosov lattice and the array of defects. 

Recently, Lu et al. [22] have reported a vortex ratchet effect for a 2D vortex system 
interacting with a substrate that has an asymmetric periodic modulation along one 
direction. Although the pinning in that system might be expected to have a ID 
character, Lu et al showed that the ratchet effect and ratchet reversals that arise are 
produced by collective vortex-vortex interactions, and also showed that the overall 
vortex lattice structure can change for drives in the opposite direction. 

The paper is organized as follows. In section 2 we review the experiment and 
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present results with Cu non-magnetic artificial defects. In section 3 we introduce the 
model used in numerical simulations. In section 4 we discuss current reversal and 
lattice configurations. Finally, we present our conclusions in section 5. 

2. Experimental method 

We have fabricated the Nb thin films on top of arrays of Cu or Ni nanotriangles grown 
on Si (100) substrates, by means of electron beam lithography, magnetron sputtering 
and ion etching techniques [iniUS]. The array aspect ratio and the triangle dimensions 
of the samples are similar to the films reported by Villegas et al. [10] : the periodicity 
of the array is 770 nm x 746 nm, the triangle base is 620 nm and height is 477 nm. 
A cross-shaped patterned bridge allows us to inject a current into the sample and to 
measure the voltage drop parallel or perpendicular to the triangle base. In all the 
experiments, the magnetic field H is applied perpendicular to the substrate. In the 
ratchet measurements an injected ac current is applied and the output dc voltage is 
recorded, whereas a commercial He cryostat is used to control the temperature. More 
fabrication and characterization details have been reported elsewhere [lOl [23] . 

The dc magnetoresistance of superconducting thin films with periodic arrays of 
pinning centers show minima when the vortex lattice matches the unit cell of the 
array [53] . These minima are sharp (strong reduction of the dissipation) and equally 
spaced (two neighbor minima are always separated by the same magnetic field value). 
Consequently, the number of vortices per array unit cell can be known by simple 
inspection of the dc magnetoresistance R{H) curves, in which the first minimum 
corresponds to one vortex per unit cell, the second minimum to two vortices per unit 
cell, and so on. Moreover, the ratio between the dimension of the pinning center 
and the superconducting coherence length governs the maximum number of vortices 
that could be pinned in each one of the pinning centers |25| . The dimensions of our 
defects allow a maximum of three pinned vortices per triangle. Therefore, increasing 
the applied magnetic field to four vortices per array unit cell leads to three pinned 
vortices and one interstitial vortex [10] . In general, we know, for selected values of 
the applied magnetic field, how many vortices are per unit cell and how many of them 
are pinned and interstitial in the ground state and for zero external force. We have 
to notice that, in this type of experiments, the periodic potentials govern the vortex 
dynamics inducing many remarkable effects, as vortex lattice reconfiguration [26] or 
vortex lattice channeling effects [57] . Finally, we have used different frequencies, up to 
lOkHz, for the input signal with identical results. This indicates that our experiments 
are done in the adiabatic regime [23l[14], i.e., the period of the signal is much longer 
than the relaxation time of the system, which is able to adopt the corresponding 
stationary state at any time. 

Figure [1] shows the experimental results of Nb films on arrays of Ni (magnetic) 
and Cu (non-magnetic) nanotriangles for n — A vortices per triangle. The results 
corresponding to Ni defects are similar to those reported in Ref. [TU] with a clear 
current reversal. The Cu non-magnetic defects exhibit a similar behavior. In this 
type of defect, pinning is weaker [28j and, consequently, rectification is lower and 
depinning forces are smaller. 
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Figure 1. Velocity vs. driving force for an applied magnetic field is H = 
1.3 X 10^^ T (n = 4 vortices per array unit cell) and temperature T/Tc = 0.99. 
The ac current frequency is 10 kHz. The sample is a Nb film 100 nm thick with 
array of 40 nm thick triangles (array periodicity 770 nm X 746 nm, triangle side 
620 nm). Dots: Nb film with array of Ni triangles. Triangles: Nb film with array 
of Cu triangles. 




3. Theoretical model 

To clarify the mechanisms behind the experimental results, we have performed 
extensive numerical simulations of vortices as a set of interacting and overdamped 
two-dimensional Brownian particles. Similar simulations have been used to study the 
DC vortex transport properties of arrays of triangular pinning antidots [29j and current 
reversal in arrays of triangular pinning defects [19], although in these references the 
effects of possible lattice reconfigurations are not analysed. 

The corresponding Langevin equation for position of vortex i reads: 

iqr.it) = - ^ V,?7™ (|r, - r, I) + Fpin„(r,) + Fext(t) +T,{t) (1) 

where 77 is the viscosity, Uvv{r) is the interaction potential between vortices, Fpinn(r) 
the force due to pinning artificial defects, Fext(i) the external force due to the apphed 
currents, and Ti{t) are white Gaussian noises accounting for thermal fluctuations: 

(r,(i) • Tjit')) = AkTr,5,j5{t - t'). (2) 

For the vortex interaction potential, we have used [15j : 

where = hc/{2e) is the quantum of flux, /ig is the magnetic permeability of vacuum, 
Kq is the zeroth-order modified Bessel function, and A is the penetration depth of the 
material, which is very sensitive to temperature for the conditions of the experiment, 
very close to the critical temperature. 

The external force Fcxt is the Lorentz force resulting from the applied current. 
As in the experiment, this force is perpendicular to the basis of the triangles and 
sinusoidal in time. 

The triangular pinning defects have been simulated using an attractive potential 
in the shape of an inverted triangular truncated pyramid, that is, null force in the 
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Figure 2. Potential acting on vortices due to pinning defects. 



middle section of the triangle, and constant force Fp pointing inward in a region of 
width d along the walls of the pinning defect (see figure [2]) ■ We have used an array of 
6x6 defects with periodic boundary conditions. 

The parameters in the simulation have been chosen similar to already reported 
values in samples of Nb films with arrays of periodic pinning centers |30j . We 
take A = 360 nm and temperature T = 8.26 K, which are typical values close to 
the superconducting critical temperature. The pinning strength Fp — 0.12 pN and 
the width of the pinning potential d = 0.7 /im have been set so that vortices start 
depinning approximately at the same applied force as observed in the experiment (see 
[To]), whereas the dimensions of the triangles in the simulation correspond to those 
in the experiment. We point out that nevertheless simulations performed with small 
variations on these parameters exhibit similar behavior. 

Finally, we have followed the suggestion made in [31] introducing a viscosity 77(F) 
depending on the strength of the external force. The reason of this dependence is 
that viscosity is in fact the result of the interaction of vortices with random intrinsic 
pinning centers scattered all over the material. To mimic the behavior of the viscosity 
measured in Ref. [31] , we have taken in our simulations the following phenomenological 
expression for the viscosity: 

r){F) = (0.005 e^iooi^(pN)l + 6 X 10"^) pN-s/m. (4) 

Notice however that the viscosity enters in the Langevin equation ([1]) as a time 
scale. Therefore, it affects the magnitude of the velocity of the vortices, but not 
the qualitative behavior of the lattice or the mechanisms that lead to current reversal. 

As mentioned above, our experimental results, in agreement with previous work 
[53] [13], show that the system is adiabatic in the region of applied frequencies, which 
allowed us to perform our simulations also in the adiabatic limit. We first obtain 
Vdc{F) curves for different values of constant positive and negative forces (dc curves). 
These dc curves are depicted in figure [3] where the absolute value of the velocity of 
the vortex lattice is plotted against the applied force. The inversion is clearly shown 
in the inset. Once the dc curves are obtained, the corresponding ac curves Vac{Fac) 
can be easily computed by performing the following numerical integration which does 
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Figure 3. Velocity (absolute value) of the lattice as a function of the applied force, 
for positive (*) and negative (o) forces. Simulation results and linear interpolation. 
Inset shows the inversion in more detail. 
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Figure 4. Average velocity of the vortices as a function of the r.m.s. amplitude 
Fac of the external force. 

not depend on the frequency of the force (along the paper Fac stands for the r.m.s. 
value of the external force) : 

Vac{Fac) = j Vdc (i^ac V2 sin(t)^ dt. (5) 

Figure |4] shows the average velocity for 144 particles and 6x6 pinning triangles 
{n = 4) and the parameters described above, as a function of the external force 
Fac- Simulations reproduce the experimentally observed current reversal, i.e., negative 
current for small forces and positive for higher forces, decaying to zero afterward. 
Furthermore, the forces applied and the corresponding velocities are both quite similar 
to those in the experiment. 

4. Discussion: vortex lattice reconfiguration and current reversal 

Consider a single particle or vortex in the triangular potential of Fig. [21 The force to 
depin the vortex when pushing in the positive y direction through the tip of the triangle 
is approximately half the force needed to depin it through the base, when pushing in 
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a) c) 

Figure 5. The ground state of the system for zero external force. 

the negative y direction. Positive or upward is therefore the "natural" direction for 
the rectification induced by triangular defects, and it is observed experimentally for 
n = 1, 2, 3 [TUj. However, as shown in Fig. [1] for n — A downward rectification is also 
observed for small forces, and the same unexpected result is obtained for n > 4 [lOj . 

What is the mechanism of this negative rectification? Previous models |10] 
consider that interstitial vortices feel an effective ratchet potential which is roughly the 
spatial inverse of the pinning potential, so they are rectified in the negative direction 
and move independently of the pinned vortices. However, a simple visual inspection of 
our simulations reveal that the lattice moves as a whole without significant shear and 
that current reversal is due to a reconfiguration of the vortex lattice (see supplementary 
multimedia material). 

In figure[ni[a), we depict the vortex lattice that minimizes the total free energy and 
has the periodicity of the rectangular array of triangles, for n = 4 and zero external 
force. We have found this configuration by simulated annealing on a cell with a single 
defect and periodic boundary conditions. Notice that the actual ground state or 
minimal energy configuration might not in principle have the periodicity of the array 
of defects, specially for weak pinning forces. However, our simulations reveal that 
this is not the case for the parameters chosen above: figure Elja) is indeed the ground 
state of the system, consisting of a slightly distorted triangular lattice (vortices form 
isosceles triangles but no equilateral) with one interstitial and three pinned vortices. 

This lattice starts to move either upward or downward if we add a force, but 
it turns out that the depinning force is stronger in the positive direction, yielding a 
negative rectification. 

The effect can be explained by considering first a rigid lattice. When pushing 
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upward (toward the positive y direction) the y component of the total force acting 
on the three pinned vortices due to the pinning is — SFpCosa (Fig. HKb)), a being 
the angle at the basis of the triangle. On the other hand, for a negative motion, 
the total force to overcome in the y direction is (2 — cosa)Fp (Fig. EKc)). In the 
case of equilateral triangles, cos a =1/2 and both depinning forces are identical. Our 
triangles are however a little bit flattened, with base 620 nm and height 477 nm. Using 
this values and Fp — 0.12 pN, the resulting depinning forces per vortex are 0.0437 pN 
(downward) and 0.0490 pN (upward). 

Although this argument is only valid for zero temperature and a rigid lattice, the 
depinning forces observed in the simulation (Fig. |4]) are very close to these theoretical 
values: we observe downward depinning for 0.0435 pN and upward for 0.0460 pN. The 
greater discrepancy for the upward motion indicates that the plasticity of the lattice 
is more relevant in the upward depinning, whereas the lattice is barely deformed in 
the downward motion. In the supplementary multimedia material one can see how 
the whole lattice remains pinned for an external force F = 0.44 pN, whereas it does 
move downward for the same negative force F = — 0.44pN with a small deformation. 
We then conclude that the negative rectification is mainly due to the shape of the 
triangular defects. 

Increasing the external force we get upward rectification. For F — 0.5 pN, the 
lattice is already depinned in the upward direction, but we observe a very interesting 
phenomenon which is the main result of the paper: the lattice undergoes a transition 
(see supplementary multimedia material) adopting the configuration depicted in 
figure [6] with only one (or two at most) pinned vortices. This new configuration 
is the result of a 90° rotation of the ground state of Fig. [5l and induces a relatively 
fast columnar motion of vortices which exceeds the motion in the downward direction 
for the same magnitude of the external force, yielding a net upward rectification. 

Therefore, our simulations reveal that the current reversal is entirely due to 
a reconfiguration of the vortex lattice. Changes in lattice configuration have been 
extensively reported in the literature, mainly in equilibrium [26l [32j [33] . In our case, 
however, the reconfiguration is a dynamical phenomenon. In spite of this, we can gain 
some insight on the problem by calculating the interaction energy of the configurations 
depicted in Figs. [5] and [6] and comparing it with the interaction energy of triangular 
and square lattices with the same density of vortices. 

Let a and b be respectively the horizontal and vertical distance between triangles 
in the array of defects. The lattice in Fig. [5] (we will call it H lattice) is given by 
n'VH 1 + rn^H 2, with n, m = 0, ±1, ±2, . . . and the following generating vectors: 

■ v„, . ( ) ; - ( «/J ) 

The corresponding generating vectors for the lattice in Fig. [SI denoted as V, are: 

vv,i = (o); ^v-,2=(^/J) (7) 

One can show that, for a square array of defects, i.e., a — b, lattice H is a 90° 
degree rotation of lattice V and the interaction energy is obviously identical. 

The ground state of the system without triangular defects is a regular triangular 
lattice T with vortex density a, whose generating vectors can be written as: 
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Figure 6. Configuration for -Fext = 0.1 pN. 



For an aspect ratio h/a = a/3/2, the vortex density is cr = 8/(a^\/3) and lattice 
H coincides with the regular triangular lattice T. For b/a — 2/-\/3, one can also prove 
that lattice V becomes regular. Therefore, any deviation (as in our experiment) from 
these aspect ratios induces a distortion in the lattices in order to conform with the 
periodicity of the array. In fact, since the aspect ratio of the array in the experiment 
and simulations b/a = 746/770 — 0.969 is closer to •\/3/2 rather than to 2/-\/3, H is 
less distorted than V. 

Another trivial way of conforming to the periodicity of the array with n = 4 
vortices per defect, is the rectangular lattice R generated by: 

a/2 \ _ / 

) ' ""^'^ - b/2 

We have evaluated the interaction energy. Eh, Ey, Et, Eji, of these four lattices 
for the values a = 770 nm, b = 746 nm, and a = 4/(o6) corresponding both to the 
experiment and simulations. The energies are very close to each other. In Fig. [7] we 
plot the energies compared to the ground state T, i.e., the ratios Eh/Et, Ey/Ex, 
Er/Et, as a function of A. As expected, the minimal interaction energy corresponds 
to the regular lattice T, followed by the less distorted triangular one, which is H for 
the aspect ratio 0.969 in the experiment. Since this aspect ratio is close to 1, the 
energies of H and V are very similar. The rectangular lattice, on the other hand, is 
the most energetic configuration and, consequently, the most unstable. Accordingly, 
we have never seen this type of configuration in our simulations. 



V_R.l = 



(9) 
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Figure 7. Interaction energy for tlie H, V, and R lattices compared to the regular 
triangular one, T: EnjEq^ (red squares), Ey/Ej- (blue triangles), En/Ej' (black 
circles). Computed numerically to convergence, including 10000 neighbours. 

We have of course to consider the pinning energy, which favors, in absence of 
external force, lattice H, becoming the ground state of the system, as observed in the 
simulations (see Fig. [5]). However, the above calculations point out that the energy 
cost, due to vortex interaction, of rotating the H configuration to V is relatively small. 
Since the pinning energy is, by construction, comparable to the interaction energy, it 
is reasonable to explain the response of the system as a competition between lattices 
H and V. 

It can also be drawn from the previous discussion that a substrate designed with 
a different geometry (for example, a lattice of defects with an aspect ratio closer to 
2 / VS) could favour the appearance of one of the lattices over the other and this might 
prevent the transition between them, affecting current reversal. 

5. Conclusions 

We have shown that current reversal in superconducting ratchets can be explained as a 
dynamical reconfiguration of the vortex lattice. We have presented numerical evidence 
of this mechanism and also showed that the interaction energy of the involved lattices 
is close enough to allow such reconfiguration. 

The present work prompts further research in several lines. First, the extension 
of experimental and theoretical studies on lattice transitions [551 [321 123] to the case 
of moving lattices. Secondly, the analysis of the interplay between rectification and 
lattice configuration. We have shown that the intensity and direction of rectification 
can be affected not only by the shape but also by the array of defects. This could help 
to design rectifiers with specific responses to external forces as well as become a tool 
to study lattice configurations. For instance, magnetoresistance, i.e., the response of 
the system to dc currents, is commonly used to reveal the underlying vortex lattice in 
superconducting films [TUl [H] • Our work indicates that the response to ac currents in 
the presence of asymmetric artificial defects can also provide valuable information on 
the lattice structure. 

Finally, similar ratchet mechanisms could be observed in other type of systems. 
The main ingredient is a lattice or complex object undergoing conformational changes 
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as a response to external forces and in order to adapt to some anisotropic substrate. 
This conformation changes wiU depend on the sign of the force and can amplify or 
modify the rectification properties of the substrate. We believe that this type of 
generic ratchet effect could be observed in systems with interacting particles, like 
colloidal suspensions or granular media. 
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